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We prove an exact solution of a multi-lane totally asymmetric simple exclusion process (TASEP) 
with heterogeneous lane-changing rates on a torus. The solution is given by a factorized form; that 
is, the TASEP in each lane can be separable because of the detailed balance conditions satisfied for 
lane-changing transitions. Using the saddle point method, the current of particles is calculated in 
a simple form in a thermodynamic limit. It is notable that the current depends only on a set of 
lane-changing parameters, not on the configuration of lanes. 


I. INTRODUCTION 


Driven diffusive systems have been studied actively in 
recent years, since they are useful for understanding var¬ 
ious phenomena in physics and biology. One of the most 
important driven particle systems, the totally asymmet¬ 
ric simple exclusion process (TASEP) [I|, was originally 
proposed as a model for describing biological transport 
phenomena, and has been applied to the modeling of 
transport processes such as vehicular traffic 0 , granular 
flow [3| , and biological transportation by motor proteins 

M. 


In some studies, the TASEPs with multiple lanes and 
lane-changing have been investigated analytically [rl-fT^ , 
but exact analyses have been performed on a few models 


In this work, we consider a multilane system with pe¬ 
riodic boundaries in two directions, and present an exact 
solution in the stationary limit. The system has K lanes 
on a cylinder, which is applicable to the problems such as 
transportation phenomena of the kinesins [l3l| along 
the 13 protofilaments placed on microtubules cylindri- 
cally [ij, [2. On the other hand, when K = 2 it corre¬ 
sponds to a simple two-lane TASEP with periodic bound¬ 
ary. Moreover, we do not limit the number of lanes, and 
thus this work will be a significant achievement for solv¬ 
ing a kind of two-dimensional exclusion process exactly. 

To construct the solution, we use the detailed bal¬ 
ance condition satisfied in lane-changing transitions in 
the model. For this characteristic, the solution has a sim¬ 
ple structure, and quantities such as density and current 
are derived in simple forms. 


II. MODEL 

We consider a two-dimensional cylindrical lattice of L x 
K sites. Each lane-f [i = 1, • • ■ ,K) has L sites (j = 
1, • • • ,L), and corresponding sites of adjacent lanes are 
connected with each other. These K lanes are arranged 
cylindrically, namely, lane-i -I- AT is identical with lane-i. 

A site-(i,j) can be either be empty (rij = 0) or occu¬ 


pied by one particle (ji^j = 1) (the hard-core exclusion). 
Ti denotes the occupation number of the j-th site in 
lane-f. The time evolution per time interval At is writ¬ 
ten as follows: 

(i) hopping 

lijOij+i —?> Oijlij+i with probability PiAt 

(ii) lane-changing 

-5> Oij, li±i,j with probability XiAt 

In each lane, a particle hops to the next site (j -T j -\-l) 
with probability PiAt{> 0) if the target site is empty (as 
in the usual TASEP). Furthermore, lane-changing tran¬ 
sitions are also defined, namely, each particle hops to the 
adjacent lane {i ^ i -\- 1 or i ^ i — 1) with probability 
XiAt{> 0) if the target site is empty as shown in FIG[T] 
Since in our model each particle is randomly updated, 
this lane-changing is permitted even if the neighboring 
site-j — 1 to the target site is occupied. 



FIG. 1. The schematic view of the system. If the neighboring 
site to its right is empty, the particle hops there j j + 1 with 
probability piAt (hopping). Moreover, if the receiving site is 
empty, it changes lanes to the neighboring lane i -|- 1 (i — 1) 
[i,j —>■ {i—l,j)] with probability XiAt (lane-changing). 


When K = 2 and each lane has only one adjacent lane 
the lane-changing is restricted to one direction. This pe¬ 
culiarity is natural and does not influence the argument 
in the following sections. 


















2 


III. EXACT SOLUTION FOR A PERIODIC 
SYSTEM 

In the following we focus on the periodic system in the 
hopping direction; namely, site-j + L is identical with 
site-j. For this system, we present an exact expression for 
the probability with which a given configuration {r^j} is 
realized, as a product form of the density weight 'I'i(Mi) 
and the configuration weight in lane-f, Here, 

Mi is the number of particles in lane-z {Mi = Tij ) and 
= Ti,!,'' ’ R.L S {0,1} is a Set of occupation 
numbers in lane-i. 

K 

i=l 

K / 1 X M, 

= ^L^N,KW_\—j ( 2 ) 

Here, each weight factor gi{{Ti^j}i) is assumed to be 1 as 
in the usual ASEP with periodic boundary [T^, and the 
density weight '^i{Mi) is defined as '^i{Mi) = {1/xi)^'- 
The normalization factor is thus written as 

K 

i=i 

( 3 ) 

by taking the sum of the weights with respect to Mi. The 
(5-function ensures that only those configurations with the 
correct total number of particles are included. 

Next, we confirm that this expression gives the exact 
solution. The exact solution for the steady state of the 
system must satisfy the master equation described be¬ 
low, and conversely, the expression satisfying the master 
equation must be the exact solution. 

0 = ^P{C) = {P{C')W{C' ^ C) - P{C)W{C C')} 

( 4 ) 

Here, C and VF(C —>■ C') indicates the configuration of 
particles and the transition probability from configura¬ 
tion C to C' respectively. 

We separate the transitions into two parts according 
to their type of motion, i.e., hopping or lane-changing. 
It is obvious that the terms for the hopping transition 
in the master equation (|31) vanish when one substitutes 
the presented solution, since each gi{{Tij}i) is the exact 
solution of the TASEP with periodic boundary in each 
lane. For transition terms of hopping in lane-i, one only 
has to consider the weight of lane-i gi{{Tij}i) since other 
terms in Eq. o are in common before and after the 
transition. 

Then we show that the rest of Eq. o, namely, 
the terms for lane-changing transitions, also vanish. 
We focus on the lane-changing transitions concerned 
with lane-i. It is sufficient to consider the neighbor¬ 
ing lanes {i — l,i,i + 1), and we write the configu¬ 
rations of these three lanes with and Mi+i 


L LK, 

Zl,n,k ^ Y ■■■ — 

Mi=0 MK=0i=l 


Mi 


L 

Mi 


particles as j}- First, let us confirm 

the correspondence of the lane-changing transitions i —> 
i + 1 and z + 1 —>• i. As shown in FIGHJ the num¬ 
ber of configurations of these transitions are the same 
where the asterisks indicate the common configurations 
among the lanes, namely, every lane-changing transi¬ 
tion {• • • , Tij = 1, • • • , Ti±ij = 0, • • • } {■■■ ,Tij = 

Oj ■ • ■ )R±i.j = I)"'"} is paired with its counterpart 

{' ’ * ( Rj — ‘ ' 5 Piljj — 1; ' ’ ' 1 ^ ' 5 — 

1; ■ • ■ ! These corresponding transitions 

balance in the master equation as explained below. Here, 
we avoid the discussion with explicit expressions expand¬ 
ing 0 because the expressions would be unnecessarily 
complicated and make us lose sight of the essence. The 
balance of lane-changing transitions concerned with lane- 
i is illustrated in FIG El The transitions not described 
in the figure are forbidden in this model (they occur 
with probability 0). We choose one arbitrary configu¬ 
ration of lanes z — l,i, and z -I- 1, 

Taking the transition between and 

broken line in the fig¬ 
ure) as an example, they cancel in the master equation 
as follows: 


P{C')W{C' C) - P{C)W{C C) (5) 



= 0 (7) 


One can understand the rest of transitions also balance 
using the same argument. Since the system is periodic, 
this balance holds for every lane and configuration. It 
should be noted that the transitions between the neigh¬ 
boring lanes satisfy the detailed balance condition. Thus, 
we have proved that the presented expression surely sat¬ 
isfies the master equation and correctly describes the sys¬ 
tem. 



* ■■■lane - i+l 


^ ■ ■ ■lane - i 


i+l 


I+1—> i 


FIG. 2. Correspondence of lane-changing transitions on a 
certain site in lanes i and z -|- 1. Asterisks indicate the same 
configuration in both transitions. 
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FIG. 3. Possible lane-changing transitions concerning lane-i for an arbitrary set of configurations Cm^ 

represents a certain configuration in lane-i, where Mi is the total number of particles in lane-i in the configuration. 


IV. THERMODYNAMIC LIMIT 


Let us discuss the thermodynamic limit (L oo) of 
the system. First, we consider the following functions. 


MM) = ( - 

\Xi 


M 


L 

M 


ZlMs)= E ••• E ]lMMi)i 

Ml—0 MjK—Oi—l 
K L 


Mi 


=n E 

i^l Mi^ 

^ / 

=n 1 


Mi 


2=1 
K 


S 

Xi 


n(-y 


= [MM 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


where MM) is the weight of lane-i and F{s) = 
I\f=i (1 + s/x^)■ Using this Zl.k{s), the partition func¬ 
tion is expressed in an integral form. 


Zl,n,k 


ds Zl,k{s) 
2m 


ds [Fjs)]^ 

2TTi 


(13) 


We evaluate the integral (fT^ in the L ^ oo limit, keeping 
N/L = p{< K) constant, by the saddle point method. 
For large L, N Eq. (IT^ is dominated by the saddle point 
of the integral denoted by s = z. Following M , we define 

(^(s) = —pins-h ln[F(s)]. (14) 

Then the saddle point is given by ((>'(z) = 0 as. 


K 


^ = E 
2 = 1 


z/Xi 

1 + z/Xi' 


(15) 


saddle point as, 

^ _ 1 _ 1 exp(L(^(z)) 

“ (27rL)i/2 |</."(z)|i/2 2 

1 1 Zl,k{z) 

(27rL)i/2 |<(."(z)|i/2 ^Af+i 


(16) 


Let us investigate the current in lane-l which is defined 
as Ji = {Tip{l — Tij+i)pi), by considering the configu¬ 
rations where Tij = l,Tij+i = 0. The weight of these 
configurations is calculated through its generating func¬ 
tion, 


K 




L / Mi 




Mi=0 

s/x 


L-2 
Mi - 1 


[U(s)]^ 


(1 + s/Xi)^ 

Therefore, by the same argument we find 

72 

Ji = t :- Pi 


Zl. 


N,K 


ds s/xi [^(^)]^ 

2772 (1 + s/Xi)^ S-^+^ 


ds iFjs)]^ 

2772 S-^ + 1 


z/Xi 


{^ + z/x^)^ 


Pi 


r.Mi 


(17) 

(18) 

(19) 

( 20 ) 


in the thermodynamic limit. Here, z is the saddle point of 
the integral again. Note that the two integrals performed 
above have common (j), and thus the saddle point is con¬ 
sistent. In a similar way, the density in lane-f pt = Mi/L 
is also calculated as, 


Pi ^ 


z/Xi 
1 + z/xi 


( 21 ) 


Moreover, the partition function is evaluated by consid- This corresponds with the expression (1151) where each pi 
ering the thermodynamic limit of Eq. (USD using this contributes to the total density. Eurthermore, it is quite 
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notable that z plays the role of a “common incoming 
rate” when we compare it with the Langmuir equilibrium 
density k/(1 + k) where k is the ratio of the attachment 
and detachment rates 0. 

The currents Ji obtained from Monte Carlo simula¬ 
tions are plotted on FIG. |4] with theoretical lines. The 
theoretical lines are obtained from Eq. dsni) after one 
finds z from Eq. (USD. It is notable that the density in 
lane i depends only on the set of lane-changing rates, and 
is independent of the configuration of the lanes. For rela¬ 
tively large leaving rate of lane-i Xi density becomes 
small, and leads to the large critical density. 

To summarize this discussion, we can also r^ard the 
dynamics as the ASEP with Langmuir kinetics Q on each 
lane with detachment rate Xi effective attachment 
rate z in the thermodynamic limit for its detailed balance 
property. 



V. SUMMARY 


In this work we have considered a cylindrical multi¬ 
lane exclusion process and presented an exact solution of 
it in the stationary limit. This solution is also applicable 
to the model with K = 2 lanes, which corresponds to a 
two-lane model often discussed in transportation prob¬ 
lems. Using the saddle point method we have derived an 
expression for the current and density, and have shown 
that simulation results with representative parameters 
well agree with the theory. As shown in EIGlU the peak 
shift has been observed, and this phenomenon might be 
seen in the actual transport process such as the biological 
transportation of motor proteins. 


The most important feature of the model is the de¬ 
tailed balance condition satisfied in the lane-changing 
transitions; and for this characteristic, the solution has 
been constructed simply and interpreted as separated 
ASEPs with Langmuir kinetics with a common rate in 
the thermodynamic limit. 


It is also intriguing that the system is solvable although 
the one-dimensional heterogeneous symmetric simple ex¬ 
clusion process (which corresponds to the dynamics in 
the i direction in this work) itself has not been solved in 
the previous works so far. Moreover, it is significant that 
an exact solution for a two-dimensional exclusion process 
has been given. 


FIG. 4. Fundamental diagram for the system with four lanes, 
L = 1000, (xi, X 2 , X 3 , X4) = (0.1, 0.2, 0.05, 0.3) and pi = p 2 = 
P3 = P4 = 0.3. Here, the hopping probabilities pi are set 
as constant to make the peaks of each curve aligned for the 
comprehensibility. The relationship of the curves depends on 
the lane-changing probability, and the relatively large Xi leads 
to the large critical density. 


In the model we have considered a system with sym¬ 
metric lane-changing rates, where decrease of particles 
in the lane corresponds to one lane-changing parameter. 
If we assume asymmetric ones, the formulation of the 
expression would be more complex, and it should be in¬ 
vestigated in future works. 
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